Pre- and early-postnatal nutrition modify gene and protein expressions of muscle energy metabolism markers and phospholipid fatty acid composition in a muscle type specific manner in sheep.

We previously reported that undernutrition in late fetal life reduced whole-body insulin sensitivity in adult sheep, irrespective of dietary exposure in early postnatal life. Skeletal muscle may play an important role in control of insulin action. We therefore studied a range of putative key muscle determinants of insulin signalling in two types of skeletal muscles (longissimus dorsi (LD) and biceps femoris (BF)) and in the cardiac muscle (ventriculus sinister cordis (VSC)) of sheep from the same experiment. Twin-bearing ewes were fed either 100% (NORM) or 50% (LOW) of their energy and protein requirements during the last trimester of gestation. From day-3 postpartum to 6-months of age (around puberty), twin offspring received a high-carbohydrate-high-fat (HCHF) or a moderate-conventional (CONV) diet, whereafter all males were slaughtered. Females were subsequently raised on a moderate diet and slaughtered at 2-years of age (young adults). The only long-term consequences of fetal undernutrition observed in adult offspring were lower expressions of the insulin responsive glucose transporter 4 (GLUT4) protein and peroxisome proliferator-activated receptor gamma, coactivator 1α (PGC1α) mRNA in BF, but increased PGC1α expression in VSC. Interestingly, the HCHF diet in early postnatal life was associated with somewhat paradoxically increased expressions in LD of a range of genes (but not proteins) related to glucose uptake, insulin signalling and fatty acid oxidation. Except for fatty acid oxidation genes, these changes persisted into adulthood. No persistent expression changes were observed in BF and VSC. The HCHF diet increased phospholipid ratios of n-6/n-3 polyunsaturated fatty acids in all muscles, even in adults fed identical diets for 1½ years. In conclusion, early postnatal, but not late gestation, nutrition had long-term consequences for a number of determinants of insulin action and metabolism in LD. Tissues other than muscle may account for reduced whole body insulin sensitivity in adult LOW sheep

Genetic architecture of gene expression in ovine skeletal muscle

In livestock populations the genetic contribution to muscling is intensively monitored in the progeny of industry sires and used as a tool in selective breeding programs. The genes and pathways conferring this genetic merit are largely undefined. Genetic variation within a population has potential, amongst other mechanisms, to alter gene expression via cis- or trans-acting mechanisms in a manner that impacts the functional activities of specific pathways that contribute to muscling traits. By integrating sire-based genetic merit information for a muscling trait with progeny-based gene expression data we directly tested the hypothesis that there is genetic structure in the gene expression program in ovine skeletal muscle. Results The genetic performance of six sires for a well defined muscling trait, longissimus lumborum muscle depth, was measured using extensive progeny testing and expressed as an Estimated Breeding Value by comparison with contemporary sires. Microarray gene expression data were obtained for longissimus lumborum samples taken from forty progeny of the six sires (4-8 progeny/sire). Initial unsupervised hierarchical clustering analysis revealed strong genetic architecture to the gene expression data, which also discriminated the sire-based Estimated Breeding Value for the trait. An integrated systems biology approach was then used to identify the major functional pathways contributing to the genetics of enhanced muscling by using both Estimated Breeding Value weighted gene co-expression network analysis and a differential gene co-expression network analysis. The modules of genes revealed by these analyses were enriched for a number of functional terms summarised as muscle sarcomere organisation and development, protein catabolism (proteosome), RNA processing, mitochondrial function and transcriptional regulation. Conclusions This study has revealed strong genetic structure in the gene expression program within ovine longissimus lumborum muscle. The balance between muscle protein synthesis, at the levels of both transcription and translation control, and protein catabolism mediated by regulated proteolysis is likely to be the primary determinant of the genetic merit for the muscling trait in this sheep population. There is also evidence that high genetic merit for muscling is associated with a fibre type shift toward fast glycolytic fibres. This study provides insight into mechanisms, presumably subject to strong artificial selection, that underpin enhanced muscling in sheep populations

Identification of co-expression gene networks, regulatory genes and pathways for obesity based on adipose tissue RNA Sequencing in a porcine model

Background: Obesity is a complex metabolic condition in strong association with various diseases, like type 2 diabetes, resulting in major public health and economic implications. Obesity is the result of environmental and genetic factors and their interactions, including genome-wide genetic interactions. Identification of co-expressed and regulatory genes in RNA extracted from relevant tissues representing lean and obese individuals provides an entry point for the identification of genes and pathways of importance to the development of obesity. The pig, an omnivorous animal, is an excellent model for human obesity, offering the possibility to study in-depth organ-level transcriptomic regulations of obesity, unfeasible in humans. Our aim was to reveal adipose tissue co-expression networks, pathways and transcriptional regulations of obesity using RNA Sequencing based systems biology approaches in a porcine model. Methods: We selected 36 animals for RNA Sequencing from a previously created F2 pig population representing three extreme groups based on their predicted genetic risks for obesity. We applied Weighted Gene Co-expression Network Analysis (WGCNA) to detect clusters of highly co-expressed genes (modules). Additionally, regulator genes were detected using Lemon-Tree algorithms. Results: WGCNA revealed five modules which were strongly correlated with at least one obesity-related phenotype (correlations ranging from -0.54 to 0.72, P <0.001). Functional annotation identified pathways enlightening the association between obesity and other diseases, like osteoporosis (osteoclast differentiation, P = 1.4E(-7)), and immune-related complications (e. g. Natural killer cell mediated cytotoxity, P = 3.8E(-5); B cell receptor signaling pathway, P = 7.2E(-5)). Lemon-Tree identified three potential regulator genes, using confident scores, for the WGCNA module which was associated with osteoclast differentiation: CCR1, MSR1 and SI1 (probability scores respectively 95.30, 62.28, and 34.58). Moreover, detection of differentially connected genes identified various genes previously identified to be associated with obesity in humans and rodents, e.g. CSF1R and MARC2. Conclusions: To our knowledge, this is the first study to apply systems biology approaches using porcine adipose tissue RNA-Sequencing data in a genetically characterized porcine model for obesity. We revealed complex networks, pathways, candidate and regulatory genes related to obesity, confirming the complexity of obesity and its association with immune-related disorders and osteoporosis

Possible mechanisms of host resistance to Haemonchus contortus infection in sheep breeds native to the Canary Islands

Haemonchus contortus appears to be the most economically important helminth parasite for small ruminant production in many regions of the world. The two sheep breeds native to the Canary Islands display distinctly different resistant phenotypes under both natural and experimental infections. Canaria Hair Breed (CHB) tends to have significantly lower worm burden and delayed and reduced egg production than the susceptible Canaria Sheep (CS). To understand molecular mechanisms underlying host resistance, we compared the abomasal mucosal transcriptome of the two breeds in response to Haemonchus infection using RNAseq technology. The transcript abundance of 711 and 50 genes were significantly impacted by infection in CHB and CS, respectively (false discovery rate <0.05) while 27 of these genes were significantly affected in both breeds. Likewise, 477 and 16 Gene Ontology (GO) terms were significantly enriched in CHB and CS, respectively (P < 1.0 × 10(−4)). A broad range of mechanisms have evolved in resistant CHB to provide protection against the parasite. Our findings suggest that readily inducible acute inflammatory responses, complement activation, accelerated cell proliferation and subsequent tissue repair, and immunity directed against parasite fecundity all contributed to the development of host resistance to parasitic infection in the resistant breed

Detection of regulator genes and eQTLs in gene networks

Genetic differences between individuals associated to quantitative phenotypic traits, including disease states, are usually found in non-coding genomic regions. These genetic variants are often also associated to differences in expression levels of nearby genes (they are "expression quantitative trait loci" or eQTLs for short) and presumably play a gene regulatory role, affecting the status of molecular networks of interacting genes, proteins and metabolites. Computational systems biology approaches to reconstruct causal gene networks from large-scale omics data have therefore become essential to understand the structure of networks controlled by eQTLs together with other regulatory genes, and to generate detailed hypotheses about the molecular mechanisms that lead from genotype to phenotype. Here we review the main analytical methods and softwares to identify eQTLs and their associated genes, to reconstruct co-expression networks and modules, to reconstruct causal Bayesian gene and module networks, and to validate predicted networks in silico.Comment: minor revision with typos corrected; review article; 24 pages, 2 figure

RNA-Seq transcriptomics and pathway analyses reveal potential regulatory genes and molecular mechanisms in high- and low-residual feed intake in Nordic dairy cattle

BACKGROUND: The selective breeding of cattle with high-feed efficiencies (FE) is an important goal of beef and dairy cattle producers. Global gene expression patterns in relevant tissues can be used to study the functions of genes that are potentially involved in regulating FE. In the present study, high-throughput RNA sequencing data of liver biopsies from 19 dairy cows were used to identify differentially expressed genes (DEGs) between high- and low-FE groups of cows (based on Residual Feed Intake or RFI). Subsequently, a profile of the pathways connecting the DEGs to FE was generated, and a list of candidate genes and biomarkers was derived for their potential inclusion in breeding programmes to improve FE. RESULTS: The bovine RNA-Seq gene expression data from the liver was analysed to identify DEGs and, subsequently, identify the molecular mechanisms, pathways and possible candidate biomarkers of feed efficiency. On average, 57 million reads (short reads or short mRNA sequences < ~200 bases) were sequenced, 52 million reads were mapped, and 24,616 known transcripts were quantified according to the bovine reference genome. A comparison of the high- and low-RFI groups revealed 70 and 19 significantly DEGs in Holstein and Jersey cows, respectively. The interaction analysis (high vs. low RFI x control vs. high concentrate diet) showed no interaction effects in the Holstein cows, while two genes showed interaction effects in the Jersey cows. The analyses showed that DEGs act through certain pathways to affect or regulate FE, including steroid hormone biosynthesis, retinol metabolism, starch and sucrose metabolism, ether lipid metabolism, arachidonic acid metabolism and drug metabolism cytochrome P450. CONCLUSION: We used RNA-Seq-based liver transcriptomic profiling of high- and low-RFI dairy cows in two breeds and identified significantly DEGs, their molecular mechanisms, their interactions with other genes and functional enrichments of different molecular pathways. The DEGs that were identified were the CYP’s and GIMAP genes for the Holstein and Jersey cows, respectively, which are related to the primary immunodeficiency pathway and play a major role in feed utilization and the metabolism of lipids, sugars and proteins. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12864-017-3622-9) contains supplementary material, which is available to authorized users

Quantitative Trait Loci Involved in Sex Determination and Body Growth in the Gilthead Sea Bream (Sparus aurata L.) through Targeted Genome Scan

Among vertebrates, teleost fish exhibit a considerably wide range of sex determination patterns that may be influenced by extrinsic parameters. However even for model fish species like the zebrafish Danio rerio the precise mechanisms involved in primary sex determination have not been studied extensively. The zebrafish, a gonochoristic species, is lacking discernible sex chromosomes and the sex of juvenile fish is difficult to determine. Sequential protandrous hermaphrodite species provide distinct determination of the gender and allow studying the sex determination process by looking at the mechanism of sex reversal. This is the first attempt to understand the genetic basis of phenotypic variation for sex determination and body weight in a sequential protandrous hermaphrodite species, the gilthead sea bream (Sparus aurata). This work demonstrates a fast and efficient strategy for Quantitative Trait Loci (QTL) detection in the gilthead sea bream, a non-model but target hermaphrodite fish species. Therefore a comparative mapping approach was performed to query syntenies against two other Perciformes, the European sea bass (Dicentrarchus labrax), a gonochoristic species and the Asian sea bass (Lates calcarifer) a protandrous hermaphrodite. In this manner two significant QTLs, one QTL affecting both body weight and sex and one QTL affecting sex, were detected on the same linkage group. The co-segregation of the two QTLs provides a genomic base to the observed genetic correlation between these two traits in sea bream as well as in other teleosts. The identification of QTLs linked to sex reversal and growth, will contribute significantly to a better understanding of the complex nature of sex determination in S. aurata where most individuals reverse to the female sex at the age of two years through development and maturation of the ovarian portion of the gonad and regression of the testicular area. [Genomic sequences reported in this manuscript have been submitted to GenBank under accession numbers HQ021443–HQ021749.

Dissecting the Genetic Architecture of Host–Pathogen Specificity

International audienc

Genetic parameters for social effects on survival in cannibalistic layers: Combining survival analysis and a linear animal model

<p>Abstract</p> <p>Background</p> <p>Mortality due to cannibalism in laying hens is a difficult trait to improve genetically, because censoring is high (animals still alive at the end of the testing period) and it may depend on both the individual itself and the behaviour of its group members, so-called associative effects (social interactions). To analyse survival data, survival analysis can be used. However, it is not possible to include associative effects in the current software for survival analysis. A solution could be to combine survival analysis and a linear animal model including associative effects. This paper presents a two-step approach (2STEP), combining survival analysis and a linear animal model including associative effects (LAM).</p> <p>Methods</p> <p>Data of three purebred White Leghorn layer lines from Institut de Sélection Animale B.V., a Hendrix Genetics company, were used in this study. For the statistical analysis, survival data on 16,780 hens kept in four-bird cages with intact beaks were used. Genetic parameters for direct and associative effects on survival time were estimated using 2STEP. Cross validation was used to compare 2STEP with LAM. LAM was applied directly to estimate genetic parameters for social effects on observed survival days.</p> <p>Results</p> <p>Using 2STEP, total heritable variance, including both direct and associative genetic effects, expressed as the proportion of phenotypic variance, ranged from 32% to 64%. These results were substantially larger than when using LAM. However, cross validation showed that 2STEP gave approximately the same survival curves and rank correlations as LAM. Furthermore, cross validation showed that selection based on both direct and associative genetic effects, using either 2STEP or LAM, gave the best prediction of survival time.</p> <p>Conclusion</p> <p>It can be concluded that 2STEP can be used to estimate genetic parameters for direct and associative effects on survival time in laying hens. Using 2STEP increased the heritable variance in survival time. Cross validation showed that social genetic effects contribute to a large difference in survival days between two extreme groups. Genetic selection targeting both direct and associative effects is expected to reduce mortality due to cannibalism in laying hens.</p